Dynamic Coulomb blockade in single— lead quantum dots 
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We investigate transient dynamic response of an Anderson impurity quantum dot to a family of 
ramp-up driving voltage applied to the single coupling lead. Transient current is calculated based 
on a hierarchical equations of motion formalism for open dissipative systems [J. Chem. Phys. 128, 
234703 (2008)]. In the nonlinear response and nonadiabatic charging regime, characteristic resonance 
features of transient response current reveal distinctly and faithfully the energetic configuration of 
the quantum dot. We also discuss and comment on both the physical and numerical aspects of the 
theoretical formalism used in this work. 
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I. INTRODUCTION 

Dynamics of open electronic systems far from equilib- 
rium has attracted considerable interest in the past few 
years, since the dynamic processes are closely related to 
and thus can be used to explore the system properties 
of interest. It is well known that electron-electron in- 
teraction plays a significant role under realistic experi- 
mental conditions [if. Therefore, understanding the elec- 
tronic dynamics in presence of electron-electron interac- 
tion is of fundamental importance to the emerging field 
of nanosciences. 

So far theoretical work has focused mostly on steady or 
quasi-steady state dynamics, which however may contain 
only partial information on the system of study. Con- 
sider, for example, a quantum dot (QD) as depicted in 
Fig-UJa), with two Zeeman levels of e-\ < ej_, in con- 
tact with a single electrode. Under the quasi-steady 
state charging condition, the applied voltage raises the 
electrode chemical potential (/z) adiabatically. The first 
charging electron occupies exclusively the energetically 
favorable spin-up level, while the second one that occu- 
pies the spin-closed double occupation state requires an 
excess energy to overcome the static on-dot Coulomb 
blockade; see Fig.QJb). Apparently, the steady-state 
study offers rather incomplete information, as the single 
occupation of spin-down state is dark here. 

In this work, we resort to transient dynamics of the 
system, by applying a time-dependent external voltage 
to the system, i.e., p varies explicitly in real time, and an- 
alyzing the resulting nonadiabatic charging process. It is 
expected that all transitions within the range of applied 
voltage will be manifested through the system dynamic 
properties. Specifically, we consider a single-lead Ander- 
son impurity QD driven by a family of ramp-up external 
voltages, by exploiting the hierarchical equations of mo- 
tion (HEOM) theory for open electronic dynamics and 
transient current Q. In Sec.|TT]we briefly describe the 




FIG. 1: (a) An interacting QD coupled to a single lead with a 
time-dependent chemical potential p(t) = /i e q+A(t). (b) The 
steady-state number of electrons on QD, N e , as a function of 
p. Other parameters adopted are (in unit of meV): ej = 0.5, 
ei = 2.5, U = 4 , T = 0.2, the QD-lead coupl ing strength 
T = 0.05 and bandwidth W = 15. See Sec. lTTTAl for details. 



theoretical and computational aspects. We demonstrate 
in Sec. Mil that the transient dynamics reveals sensitively 
and faithfully the QD energetic configuration. Section 
IIVI concludes this work. 



II. METHODOLOGY 

We implement a HEOM formalism, which is for- 
mally exact for dynamics of an arbitrary non-Markovian 
dissipative system interacting with surrounding baths 
0, S BL H- From the perspective of quantum dissipa- 
tion theory, the electronic dynamics of an open system is 
mainly characterized by the reduced system density ma- 
trix p(t) = tr B p T (t) and the transient current through 
the coupling leads. Here p T (t) is the total density ma- 
trix of the entire system, and tr B denotes a trace over all 
lead degrees of freedom. The final HEOM is cast into a 
compact form of 



p n = -[iC + j n (t)]p n +p i n } +p l n +} . 



(1) 
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Here, C ■ = [H (t), ■ ] is system Liouvillian (setting H = 1). 
The basic variables of Eq. |T]) are p(t) and associated aux- 
iliary density operators (ADOs) p n {t), where n is an index 
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set covering all accessible derivatives of the Feynman- 
Vernon influence functional [6|. 

For a single-level QD coupled to a single lead, the lead 
correlation function can be expanded by an exponential 
series via a spectral decomposition scheme @, B Q ■ In 
this case n involves an n-fold combination of (a, s, m), 
which characterizes the exponential series expansion with 
a = ±, s the spin, and m the index of exponents, respec- 
tively. Therefore, in Eq. fll} p n \n =0 = p(t) with 7 n |fi = o = 
Pn~ } U=o = 0; p n |n>o is an ADO at the n th -tier; 7 n (£) col- 
lects all the related exponents along with external volt- 
ages, to p n (t); and p„~ } and p„ +i are the nearest lower- 
and higher-tier counterparts of p n , respectively. In par- 
ticular, the l st -tier ADOs, p n (t)|fi=i = p1 m {t), determine 
exclusively the transient current of spin-s through the 
lead as 



r ciy 

[2] 



Is(t) 



N s p T (t) 



-2Im^tr[a s p+ m W] • (2) 



N s is the lead total-occupation operator of spin s; a s 
is the system annihilation operator of spin s; and tr T 
and tr denote trace operations over the total system and 
reduced system degrees of freedom, respectively. Since 
there is only one coupling lead, the displacement cur- 
rent at every time t is sim ply — y\ I s (t); thus the Kirch- 
hoff 's current law retains [lCj, lll| . The present HEOM- 
bascd quantum transport theory admits classical geomet- 
ric capacitors or capacitive coupling to a gate electrode. 
The associated displacement current can then be evalu- 
ated readily, leading to the conservation of total current 
[ToL [Til [l2| . Gauge invariance [1(1 [HI G3 is also guaran- 
teed at all time, as it can be seen from the form of the 
th -tier of the HEOM @. 

It has been proved Q that for a noninteracting QD, the 
HEOM (fTJ) is exact with a finite terminal tier of n max = 2; 
while for an interacting QD, in principle an infinite hi- 
erarchy is required for the exact transient dynamics. In 
practice a truncation scheme is inevitable. In this work a 
straightforward truncation scheme is adopted; that is to 
set all the higher-tier ADOs zero: Pn|n>Af tru „ — 0, with 
A^ti-un the preset truncation tier. This scheme leads to a 
systematic improvement of results as A trun increases, i.e., 
by inclusion of higher-tier ADOs, as verified by extensive 
numerical tests. For a moderate system-lead coupling 
strength T < 5T, quantitatively accurate (if not exact) 
dynamics is obtained with A trun = 2, where T is the tem- 
perature in the unit of Boltzmann constant (i.e., Ub = 1)- 
A higher A trun means a drastic increase in computational 
cost. Therefore, considering the tradeoff between numer- 
ical accuracy and computational efforts, all the following 
calculations are carried out with A trun ~ 2 by confining 
the ratio T/T less than 5; see Sec. II VI for further discus- 
sions. 

The composite system of interest is described by the 
Anderson impurity model [LH ]. with its Hamiltonian H T 
expressed as 



H represents the interacting QD of our primary interest: 



H = e T n T +e i n i + U h^h i . 



(4) 



Here, e s is the energy of the spin-s (f or J.) state, 
aja s the system occupation-number operator of 



,„ s = a| 

spin s, and U the on-dot Coulomb interaction strength. 
In Eq. h B = J2k e ks h ks describes the nonin- 
teracting coupling lead, with tk a being the energy of 
its single-electron state k of spin s, and hks = d\ s dk s 
the lead occupation-number operator. The QD-lead 
coupling is given by H SB = J2k d ls a s + H - c -' 

where tk s is the coupling matrix element between the 
lead state k and the QD-level, both with spin s. The 
widely used Drude model is adopted to capture the cou- 
pling lead, i.e., a Lorentzian spectral density function of 
J(e) = TW 2 /(e 2 + W 2 ), with W characterizing the lead 
bandwidth. 

We consider explicitly the scenario depicted in 
Fig.rjTa). A ramp-up voltage V(t) is applied to the cou- 
pling lead at t = 0, which excites the QD out of equilib- 
rium. The lead energy levels are shifted due to the volt- 
age, A(t) = — eV(t), with e being the elementary charge, 
and so is the lead chemical potential p(t) = p cq + A(t). 
/i eq is the equilibrium lead Fermi energy, which is set to 
zero hereafter, i.e., p cq = 0. Under a ramp-up voltage, 
A(t) varies linearly with time until t = r, and afterwards 
is kept at a constant amplitude A: 



A(t) 



At/r, 
A, 



0<t<T 
t > T 



(5) 



H„ — H + Ht> 



(3) 



By tuning the duration parameter r, the family of ramp- 
up voltage covers three distinct regimes: (A) The adia- 
batic limit at r — > oo; (B) the intermediate range with a 
finite r; and (C) the instantaneous switch-on limit cor- 
responding to r — » + . These three regimes are indi- 
vidually explored and elaborated in the subsections of 
Sec. [mi 

The work flow of our numerical procedures are briefly 
stated here, while details to be published elsewhere, (a) 
The lead correlation function is expanded by exponential 
functions, and hence establishes Eq. (fT]). (b) The equi- 
librium reduced density matrix and ADOs in absence of 
external voltages, {p„ q }, are obtained by setting both 
sides of Eq. (fTJ to zero. The linear sparse problem is 
solved by the biconjugate gradient method [15|. (c) The 
real-time evolution driven by V(t) is solved via a direct 
integration of E q. (fTJ) by employing either the Chebyshev 
propa gato r [la , [l7| or the 4 th -order Runge-Kutta algo- 
rithm [15| . The outcomes of step (b) are adopted as the 
initial conditions for the time evolution in step (c), i.e., 
p n (t = 0) = p° q . For every time t, the transient current 
I sit) is evaluated via Eq. It normally takes about 1 
days for a typical time evolution to 100 ps with a time 
step of 0.02 ps with a single Core 2 Duo processor. 
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III. TRANSIENT RESPONSE CURRENT AND 
ITS RELATION TO THE SYSTEM ENERGETIC 
CONFIGURATION 

A. The adiabatic limit 

As r approaches oo, the infinitely slow application of 
external voltage results in quasi-static electronic dynam- 
ics for the open QD system. Consequently, at every time 
t the QD can be deemed as in the quasi-equilibrium de- 
termined by fi. This adiabatic charging scenario is intro- 
duced in Sec.[H which provides only partial information 
on the QD energetic configuration. 

The QD physical subspace is completely spanned by 
the following four Fock states, i.e., |0) (vacancy), |f) 
(single spin-up occupation), |j) (single spin-down occu- 
pation), and \d) (spin-closed double occupation). The 
QD gains an energy of e s as it transits from |0) to |s), 
or e s + U from \s) to \d), respectively. Here s labels 
the spin direction opposite to s. Each value of \x de- 
fines a steady state of zero current for the QD, together 
with the number of electrons, N e , residing on it. Shown 
in Fig.QJb) is N e as a function of fi for a specific case 
of /U eq < e-f < ej_. The two plateaus corresponding to 
quantized N e are separated by the magnitude of U, and 
the vertical step prior to each plateau represents a Fock- 
state transition which adds one more electron to the QD. 
The two steps at e-f and ej_ + U correspond to the tran- 
sitions |0) <-> |t) and |t) «-» \d), respectively; while the 
other two transitions, |0) <-> \i) and | j) «-> \d) associated 
with the excitation energies e^ and ef + U, are missing 
from Fig.QJb). This is because /i is considered to increase 
adiabatically; therefore, the first electron occupying the 
QD exclusively enters spin-up level that is energetically 
more favorable, and consequently, the second electron 
occupation requires an excess energy of U to overcome 
the static Coulomb blockade. We note that steady states 
offer rather incomplete information on the system of in- 
terest. To acquire further knowledge, it is inevitable to 
go beyond the adiabatic charging limit and explore the 
transient regime. 



B. The finite switch— on duration regime 

We now turn to cases of a finite switch-on duration 
< t < oo, in which the lead level shift A(t) is ex- 
emplified in the inset of Fig.^b). The response current 
emerges at t = immediately after the QD is excited out 
of equilibrium, and then fluctuates as the lead chemical 
potential \x sweeps over the various QD Fock states. After 
t = r, the current vanishes gradually as the reduced sys- 
tem is drawn towards the steady state with fj,(t > r) = A. 
In Fig.rSJa) we plot the transient currents for various val- 
ues of r. Peaks associated with different resonances are 
observed directly in the time domain, at its linear map- 
ping to A(t) during the period of t < r. Especially, the 
calculated I(t) for r = 50 ps resolves distinctly four peaks 
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FIG. 2: (a) Transient currents driven by a ramp-up voltage. 
The lines represent different values of turn-on duration t, 
and are separated at t — by 3nA for clarity, (b) Transient 
current of each spin direction for r = 50 ps. The ramp-up 
voltage is depicted in the inset. Other parameters are (in 
unit of meV): e T = 0.5, e L = 2.5, U = 4, V = 0.05, W = 15, 
T = 0.2, and A = 10. 




FIG. 3: (a) Spin-up and (b) spin-down transient currents in 
response to a ramp-up voltage for different U. Same param- 
eters are adopted as in Fig.[2jb). 



labeled by numbers from 1 to 4, which are responsible 
for all the four resonances consecutively activated by the 
ramp-up voltage with an increasing excitation energy of 
e-f, ej,, e-f+f/, and ei+U, respectively. Figure[2jb) depicts 
I s (t) for both spins with r = 50 ps. Within the period of 
t < t, the major peaks in I\(t) (peak-1) and Jj,(£) (peak- 
4) start to form as the time-dependent fi(t) matches e-f 
and e± + U, respectively; and they remain discernible even 
for a rather large r; see the line for r = 90 ps in Fig.[^a) . 
However, the amplitudes of the two minor peaks, peak- 
3 in 7-f (i) and peak-2 in I± (t) , diminish rapidly as the 
ramp-up duration r is lengthened. In the adiabatic limit 
(r —> oo), the minor peaks become completely invisible, 
since the spin-up and spin-down QD-lcvels are charged 
sequentially due to the Coulomb blockade; see Fig.QJb). 

Influence of Coulomb interaction on the transient cur- 
rent is illustrated by Fig. [31 Cases of different U ranging 
from a noninteracting scenario to a strongly blockaded 
QD are studied. For [7 = 0, there is only one peak 
in I s {t) for either spin, since the excitation energies e s 
and e s + U are indistinguishable. The peak is split into 
two with an increasing displacement as U is augmented. 
For either spin direction, the first peak always sticks to 
its original position as in U = 0, while the second one 
emerges later in time with a larger U. For U as large as 
5 meV, the two peaks in If (t) are almost completely sep- 
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FIG. 4: Transient current for (a) spin-up and (b) spin-down 
electrons in response to a ramp-up voltage with A = lOmeV, 
U = 4meV, and r = 50 ps. The lines represent different 
temperatures, and separated at t — by 5nA for clarity. 
Other parameters are same as in Fig.[2{b). 



arated. This confirms our attribution of the peaks to the 
various resonances: the first peaks in if (t) (peak-1) and 
1 1 (t) (peak-2) correspond to the resonances fi = ej and 
(i — e|, respectively; and the second ones in if (i) (peak- 
3) and if (t) (peak-4) are associated with [i = ef + U and 
fi = e| + U , respectively. 

Thermal effect is also explored. I s (t) under different 
lead temperatures are plotted in Fig. [4] A series of satel- 
lite oscillations appears following each resonance peak in 
I s {t) at a low T, which is resulted from a time-varying 
phase factor, exp[i J T dt A(t)], of the nonequilibrium lead 
correlation function. The oscillation frequency is thus 
determined by the magnitude of A(t). For a ramp-up 
voltage, A(t) increases monotonically during < t < t. 
It is thus reasonable to find the frequency of satellite os- 
cillations grows with time, and the decaying amplitude 
of the oscillations is due to the dissipative interactions 
between the QD and the semi-infinite lead. 

To summarize, with a finite switch-on duration r, the 
QD energetic configuration is resolved distinctly (some- 
times completely) through the resonance peaks of re- 
sponse current in real time, in obvious contrast to the 
adiabatic charging limit. 

C. The instantaneous switch— on limit 

As r approaches zero, a ramp-up voltage becomes 
asymptotically a step pulse, i.e., A(t) — A0(i), with 
0(i) being the Heaviside step function. Upon the instan- 
taneous switch-on of external voltage, both the spin-up 
and down QD-levels are activated simultaneously, pro- 
vided the voltage amplitude A is sufficiently large. It is 
intriguing to see in this limit how the response current 
would be related to the QD energetic configuration. 

In Fig. [5] (a) and (b) we plot the calculated if (t) and 
Il(t) in response to a step voltage of A = lOmeV, re- 
spectively. Three QDs of different ef are studied. The 
energetics scenario for the three QDs are /x eq = ef < ej_, 
/i oq < e-f < ex, and ef < /i eq < ej, respectively. Note that 
for all cases are in a large applied voltage regime since 




time (ps) time (ps) 

FIG. 5: (a) if(i) with ef = and (b) with three values 

of e-f in unit of meV. Other parameters are (in unit of meV): 
e x = 2.5, U = 4.5, W = 5, T = 0.05, T = 0.01, and A = 10. 



A > e s and A > e s + U, which ensure all the four QD 
Fock states accessible energetically. For both spin direc- 
tions, the transient current shows an initial overshoot- 
ing, and then oscillates on top of a decay as the reduced 
system evolves into a steady state with the lead chem- 
ical potential fi = A. These rapid oscillations are not 
artifacts. Their presence reflects the electronic dynam- 
ics of QD driven by an external voltage of large ampli- 
tude. Upon a small perturbation on ef of 0.1 meV, if (t) 
scarcely changes as those with ef = ±0.1meV almost 
overlap with the e-f = case shown in Fig.[5ja); while 
Il(t) displays drastic deviations among the three cases, 
in terms of the overshooting amplitude, the oscillation 
frequency, as well as the decaying rate. 

To understand the physical origin of the above obser- 
vations, frequency-dependent current spectrum is calcu- 
lated, i.e., I s (uj) = ^[Isit)] with T denoting the conven- 
tional Fourier transform. The results corresponding to 
the three QDs are depicted in Fig. [51 For all three cases, 
two characteristic frequencies are revealed in Re [if (a;))], 
i.e., a major kink at A — ef and a minor wiggle around 
A — ef — U. The line shape of Re [if (a;)] around the 
frequency wt = A — ef much resembles that of a single- 
level QD [ill. For ef < /i cq , Re [if (a;)] exhibits a peak 
at cjf ; while for ef > fi oq , a dip shows up. In the case of 
ef = Mcqi R- e [A ( a; )] around Wf is described by the func- 
tion form of ln(a; 2 + l)/a;, with x = 2(u>— u>i)/T. However, 
the profiles of Re[if(w)] are significantly different among 
the three cases. To be specific, for ef = /i eq , two promi- 
nent dips show up at frequencies A — ex and A — ex — U, 
respectively; for ef > /i oq , the dip at A — ex — U is consid- 
erably suppressed, which is associated with the transition 
|f) <-* \d); and for ef < fj, oq , it is the dip at A — e that 
is almost completely quenched, which indicates an inac- 
tive transition |0) «-> ||). Hence it is demonstrated the 
response current spectrum varies sensitively to the QD 
energetic configuration. 

The presence (or absence) of a resonance signature 
in Re[i s (ct>)], as well as the issue of peak versus dip at 
cjf, are closely related to the equilibrium occupancy of 
the spin-up state prior to the switch-on of the voltage. 
Specifically, in the case of ef = 0.1meV(> /i oq ), the 
initial population of spin-up electrons on the QD is as 
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FIG. 6: Re[I T (w)] (left panels) and Re[7j,(o;)] (right panels) 
for three specified values of e-j (unit of meV) in the top, mid- 
dle, and bottom panels, respectively, in relation to /i cq . Other 
parameters are those adopted in Fig.[S] and same for all pan- 
els. Insets in (a) and (b) magnify the peaks at 5.5 meV and 
lOmeV, respectively. 



small as 0.08, and the subsequent incoming spin-down 
electrons driven by the step voltage would feel scarcely 
any Coulomb repulsion due to the almost vacant spin-up 
level. This thus leads to the prominent dip at A — and 
the barely recognizable one at A — ej, — U in Re[/|(w)]; see 
Fig.EJb). Whereas in the case of = —0.1 mcV (< /i cq ), 
the spin-up level is nearly fully occupied at t — 0, i.e., 
_/Vf(0) = 0.92. Therefore, the majority of spin-down elec- 
trons injected afterwards need to acquire an excess en- 
ergy of U to overcome the Coulomb blockade, and this is 
why the dip at A — ej, becomes trivial; see Fig.^f). For 
£| = ^ cq , the spin-up level is half occupied (TVj = 0.5) 
prior to the switch-on of voltage. In such a circumstance, 
both the resonant transitions |0) <-> |J,) and |T) «-> \d) con- 
tribute equally to the response current, and hence give 
rise to the two dips of similar depth at frequencies A — e_[ 
andA — ej.— U; see Fig.[6]Jd) . With e-f >0.1meV, the res- 
onance signature in Re[/-|-(w)] at A — e-f becomes a pure 
dip, and the minor dip in Re[J|(w)] at A — — U vanishes 
completely; while with £| < — O.lmeV, Re[7|(w)] shows 
a pure peak at A — ej, and the minor dip in Re[if (a;)] at 
A — becomes invisible. 

This thus confirms that in the instantaneous switch- 
on limit, the w-dependent response current spectrum re- 
solves sensitively and faithfully the QD energetic config- 
uration through the characteristic resonance signatures. 
To our surprise, Re [/j.(w)] resolves a third peak at u>-\\ see 
the tiny bump in Fig.^b) and its inset. It possibly arises 
due to the higher-order response, and the mechanism of 
its presence requires further understanding on the tran- 



sient dynamic processes. At a lower temperature, all the 
resonance signatures in the current spectrum are accen- 
tuated (not shown), i.e., the peaks (dips) become higher 
(deeper) and sharper. In the real time, this amounts to 
an enhanced oscillation amplitude of transient current. 



IV. CONCLUDING REMARKS 

To conclude, we have investigated the transient elec- 
tronic dynamics of a single-lead interacting QD driven 
by a family of ramp-up voltages. We have found that 
(a) in the adiabatic charging limit, the quasi-steady dy- 
namics provide only incomplete information on the QD 
energetic configuration; (b) With an appropriately cho- 
sen switch-on duration, all transitions among QD Fock 
states can be resolved distinctly through the resonance 
peaks of transient current in real time, and hence pro- 
vide comprehensive (sometimes complete) information on 
the QD energetic configuration; (c) In the instantaneous 
switch-on limit, the frequency-dependent responses cur- 
rent spectrum reveals sensitively and faithfully the QD 
energetic configuration through the characteristic reso- 
nance signatures. This work thus highlights the signifi- 
cance and versatility of transient dynamics calculations. 

In the sequential tunneling regime where T <^ T (e.g., 
Figs. [IH3|), the present HEOM with truncation tier of 
-^trun = 1, which amounts to the conventional quan- 
tum master equation methods with memory will be 
sufficient. When r ~ T, calculations with iVtrun = 1 
would produce qualitatively correct static and dynamic 
Coulomb blockade for a single-lead QD, because the 
basic physics does not involve the cotunneling mecha- 
nism. However, to obtain quantitatively accurate results, 
-^trun — 2 is needed. Extensive numerical tests indicate 
that Antrim = 2 can support up to V ~ 5T (e.g., Figs. [5] 
and [5]). It has been shown [2| that the present HEOM 
theory with N tTun = 2, which properly treats the cotun- 
neling dynamics, is equivalent to a real-time diagram- 
matic formalism (ljjj : thus the Kondo problems can be 
addressed. However, a higher truncation tier is numeri- 
cally needed for convergency when T<T. The quantita- 
tively accurate results presented in this work (for T < 5T) 
are useful for further development of efficient but maybe 
approximate computational schemes. 

Investigating and/or manipulating transient electronic 
dynamics provide(s) sometimes a unique way to achieve 
physical properties and novel functions of mesoscopic 
many-body systems. However, exact theoretical results 
for transient quantum transport are rare. Weiss et al (20| 
have recently proposed an iterative real-time path inte- 
gral approach, which is numerically exact but expensive 
especially when the applied bias voltage is beyond some 
simple forms such as a step function. The present HEOM 
theory is formally exact and provides an alternative ap- 
proach. It is applicable to arbitrary time-dependent ex- 
ternal fields applied to leads and/or to the system, with- 
out extra computational cost. However, computational 
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effort increases dramatically if the converged results re- 
quire a high truncation tier of the hierarchy Q • 
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